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Abstract 

This  report  documents  two  novel  improvements  in  the  fimte-dimensional  nonlinear  dynamical  rnod- 
pling  of  marine  thrusters.  Previously  reported  models,  which  fail  to  capture  many  of  the  characteristic 
nonlinear  reponses  that  occur  during  unsteady  operations,  assume  that  the  lift  and  drag  forces  on  the 
propeller  blades  are  proportional  to  the  sine  and  cosine  of  the  angle  of  attack  where  the  angle  of  attack 
is  a  function  of  the  axial  flow  velocity  and  the  propeller’s  angular  velocity. 

We  have  found  that  the  lift  and  drag  forces  are  not  sinusoidal.  We  have  also  incorporated  the  effects 
of  rotational  fluid  velocity  and  inertia  on  thruster  response.  The  force  curves  and  model  parameters  are 
identified  using  experimental  data  from  the  load  cell  and  acoustic  doppler  current  meters.  The  accuracy 
of  the  model  is  determined  by  comparing  experimental  performance  with  numerical  simulations.  The 
results  indicate  that  thruster  models  with  nonsinusoidal  lift  and  drag  curves  provide  superior  accuracy 
in  both  transient  and  steady-state  response.  Incorporating  rotational  fluid  velocities  into  our  model 
gave  an  insignificant  improvement  for  our  case.  However,  rotational  fluid  flow  may  be  important  for 
other  types  of  thrusters.  The  research  performed  under  this  grant  was  reported  in  [9,  4, 14,  3,  5]  and  is 
referenced  at  the  end  of  the  text. 


1  Introduction 

Recent  advances  in  underwater  position  and  velocity  sensing  enable  real-time  cm-precision  position  measurements 
of  underwater  vehicles  [19,  6,  11,  7,  2,  17].  With  these  advances  in  position  sensing,  our  ability  to  precisely  control 
the  hovering  and  low-speed  trajectory  of  an  underwater  vehicle  is  limited  principally  by  our  understanding  of  ( «)  the 
vehicle’s  dynamics  and  («)  the  dynamics  of  the  bladed  thrusters  commonly  used  to  actuate  dynamically-positioned 
marine  vehicles.  This  paper  addresses  the  latter  problem.  Recent  results  indicate  that  the  transient  (unsteady) 
dynamics  of  marine  thrusters  can  be  approximated  by  a  simple  nonlinear  finite-dimensional  lumped-parameter 
dynamical  system  [18,  1,  12,  8,  16].  In  [8]  the  authors  report  a  nonlinear  thruster  dynamics  model  based  on  the 
motor  electro-mechanical  dynamics  and  thin-foil  propeller  hydrodynamics.  In  this  model,  the  propeller  and  fluid 
dynamics  are  approximated  by  a  two-dimensional  second  order  nonlinear  dynamical  system  with|state  variables  of  (t) 
axial  fluid  velocity  and  (ii)  propeller  rotational  velocity.  We  will  refer  to  this  model  as  the  “axi^  flow  model”.  In  [16] 
the  authors  report  experiments  that  corroborate  the  utility  of  the  axial  flow  model,  but  also  identify  discrepancies 
between  experimental  thruster  transient  response  and  that  predicted  by  the  model. 

This  paper  is  organized  as  follows:  First  we  review  a  previously  reported  thruster  model.  Second,  we  extend  this 
model  to  incorporate  the  effects  of  rotational  fluid  velocity  and  inertia  on  thruster  response.  Third,  we  report  a 
novel  method  for  experimentally  determining  non-sinusoidal  lift/drag  curves.  The  models  are  evaluated  by  comparing 
experimental  data  with  numerical  model  simulations.  The  data  indicates  that  thriKter  models  incorporating  both 
enhancements  provide  superior  accuracy  in  both  transient  and  steady-state  response. 

♦nrAgflT'V.aiigli  is  witk  Deep  Submergence  Lab  at  Woods  Hole  Oceanographic  Institution,  Woods  Hole,  M A,  02543,  USA, 
email:  mgrosenbaugh@whoi.edu.  Whitcomb  is  with  the  Department  of  Mechanical  Engineering,  Johns  Hopkins  University, 
Baltimore,  Maryland,  21218,  USA,  email:  llw@jhu.edu. 
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Motor  Current 

Flow  Velocity  Vector 

Axial  Flow  Velocity 

Rotational  Flow  Vel. 

Rotational  Flow  Vel.  Upstream 
Rotational  Flow  Vel.  Downstream 
Pressure 

Prop.  Rotational  Vel. 

Lift  Force 

Drag  Force 

Volumetric  Body  Force  Vector 
Hydrodynamic  Force 

Hydrodynamic  Torque 

Thrust 

Torque 

Angle  of  Attack 

Angle  of  Incidence 

Table  1:  Nomenclature 


2  Experimental  Setup 

The  data  were  obtained  with  a  newly  constructed  test  facility  capable  of  high  bandwidth  measurement  of  6-DOF 
thruster  forces  and  torques,  propeller  position,  and  3-DOF  fluid  velocity^.  The  forces  and  torques  were  all  measmed 
and  logged  at  lOOOHz.  The  force  and  torque  data  was  low-pass  filtered  with  5^^  order  zero-phase  acausal  filter 
with  cutoff  frequency  of  25Hz  to  suppress  artifacts  of  the  test  stand’s  50Hz  fundamental  vibration  mode.  All  data 
(except  force  and  torque)  are  presented  unflltered. 

3  Thruster  Dynamics  With  Axial  and  Rotational  Flow  Model 

This  section  first  reviews  a  previously  reported  finite  dimensional  thruster  model  which  considers  only  axial  fiuid 
flow  [8],  then  presents  a  more  general  model  which  includes  the  effects  of  both  rotational  and  axial  flow.  Both  models 
are  based  on  the  assumption  of  inviscid  and  incompressible  flow  with  no  velocity  component  in  the  radial  direction. 
Using  the  continuity  equation  (1)  and  the  standard  Euler  equations  (2)  [15,  13],  we  have 

VU  =  0  .  (1) 

^  +  {V*V)V  =  --Vp  +  -F.  (2) 

dt  .  ^  Q  Q  . 

The  terms  of  (1)  and  (2)  are  defined  in  Table  1.  In  the  equations  we  have  omitted  explicit  ^time  dependence  of 
variables  as  defined  in  Table  1. 

We  model  the  propeller  as  an  infinitely  thin  disc  in  the  the  center  of  the  ducted  thruster  that  applies  a  force  and  a 
torque  to  the  fluid  flow. 

F  =  0  •  Cr  +  •  Q  •  6{z)  -ee  +  Ao-F-  6{z)  -Cz  (3) 

The  reference  frame.  Figure  3,  is  in  cylindrical  coordinates,  and  S(z)  denotes  the  standard  dirac  impulse  function. 
Using  (3),  (2)  can  be  simplified  to  the  following  three  equations  in  integral  form: 

0  =  0 

‘For  a  detailed  description  of  this  test  facility  the  reader  is  referred  to  [4]  and  http://robotics.me.jhu.edu. 
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Figure  1:  Diagram  of  coordinate  frame  and  placement  of  the  disc  actuator  in  the  duct. 
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Then  applying  the  steady  Bernoulli  equation  and  linear  momentum  theory  [18,  8,  16]  to  large  control  volumes  up 
and  downstream  the  duct,  (4)  can  be  rewritten  as 


dVz 

'dt' 
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(5) 


3.1  Motor  Model 

Our  Thrusters  employ  a  direct  drive  DC-Brushless  Motor  driven  by  a  pulse  width  modulated  (PWM)  power  amplifier 
operating  in  current  control  mode.  The  motor  mechanical  dynamics  can  be  modeled  as 

Irnech  0  ’  im  Qload  fiuCtion(Dp,'op) 


where  Qioad  denotes  the  motor  shaft  load. 

3.1.1  Experimental  Determination  of  Motor  Parameters 

The  motor  has  four  parameters  Imech,  h,  fc/i  and  kfo-  The  torque  constant  h  was  provided  by  the  winding 
manufacturer.  A  least-square  value  for  the  mech£mical  inertia  of  the  system,  Imech,  was  computed  from  in-air 
experimental  data.  Careful  identification  of  the  plant  showed  that  the  friction  terms  are  dom^ated  by  linear  and 
static  friction.  Figure  2  shows  steady  state  propeller  velocity  versus  commanded  torque  in  aif  with  no  shaft  load 
applied.  Table  2  shows  the  resulting  motor  parameter  values.  In  order  to  avoid  numerical  problems  in  the  subsequent 
simulations,  static  friction  is  modeled  with  atan()  instead  of  the  discontinuous  sgn().  The  following  friction  model 
is  used  without  exception  throughout  the  paper:, 

2 

friction(f2prop)  =  *  ^prop  +  fc/0  *  —atan{20  •  flprop)  (7) 

linear  friction  static  finction 
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Figure  2:  Steady  State  shaft  velocity  versus  torque  in  air.  Motor  friction  is  dominated  by  linear  and  stick 
friction 

3.2  Propeller  Model 

The  following  lift  and  drag  functions,  [8],  are  derived  from  thin  airfoil  theory; 

D  =  \l2-Q-M-\V\^-fDm 

The  direction  of  lift  and  drag  is  perpendicular  and  parallel,  respectively,  to  the  direction  of  the  incident  flow.  Both 
are  related  to  axial  force  and  motor  torque  load  through  a  scaled  rotation  matrix. 


F' 

0  \ 

/  —  sin  a  .  —  cos  a  \ 

L 

Q. 

~\0 

0.7it  J 

cos  a  —sin  ay 

D 

3.3  Axial  Flow  Model 


In  [8]  the  authors  propose  a  single  term 
flow 


Fourier  approximation  for  lift  and  drag  functions  considering  only  axial  fluid 

/iO)  ~  Cxmoa:  '  sin(2^)  ClO^ 

fDW)=CDma^-{l-COSm)  ^  ’ 


where  the  angle  of  attack  ^  is  computed  as 


Q!  —  OiiCLTt2{^  0.7.R  Vjj) 

,0  =  $-a-  f. 


(11) 


From  (5)  the  axial  fluid  velocity  varies  as 


Vz 


\Vz\-Vz- 


(12) 


To  experimentally  identify  parameters  for  this  model  we  ran  combined  bi-directional  step  tests  in  which  the  torque 
command  to  the  motor  controller  was  instantly  reversed.  The  analjftical  model,  (6),  predicts  that  when  im  =  0  and 
Qioad  =  0  then  flprop  will  converge  to  zero  in  ^te  time,  with  the  rate  of  convergence  determined  by  the  parameter 
/.  This  anticipated  effect  is  clearly  observed  in  the  actual  experiments,  for  example  in  Figure  3  at  t  =  11s.  Figure 
3  shows  the  propeller  rotational  velocity  of  a  2.15iNrm  torque  step  and  its  reversal.  The  figure  compares  simulations 
with  different  axial  flow  length  values  to  experimental  results.  For  I  =  0.12m  the  zero  velocity  time  for  simulation 
and  experiment  match. 


Model  Independent  Parameters 


Parameter 

Value 

Description 

How  Obtained 

1 

0.12m 

duct  length,  axial  flow  length 

experiment 

Iwl 

0.12m 

outer  rotational  flow  length 

simulation 

lw2 

0.24m 

inner  rotational  flow  length 

simulation 

Ixv 

0.30m 

rotational  flow  length 

simulation 

R 

0.123m 

propeller  radius 

measured 

Ao 

0.0475m=* 

duct  cross  section  area 

measured 

# 

0.51rad 

pitch  angle 

measured 

Q 

1000  ft 

density 

property 

kt 

0.5392^ 

motor  torque  constant 

manufacturer 

0.005913^^ 

linear  friction  coefficient 

experiment  least  squares  fit 

*/o 

0.07278iVm 

stick  friction 

least  squares  fit 

7mech 

0.0014804^^ 

motor, shaft  and  prop,  inertia 

experiment  least  squares  fit 

Model  Dependent  Parameters 

Axial  Flow  model 

Ki 

axial  flow  form  factor 

experiment  least  squares  fit 

C/Lmax 

max.  Lift  coefficient 

experiment  least  squares  fit 

CJDmax 

max.  drag  coefficient 

experiment  least  squares  fit 

\  Rotational  Flow  model  | 

Ki 

axial  flow  form  factor 

experiment  least  squares  fit 

K2 

1.1415. 

rotational  flow  form  factor 

experiment  least  squares  fit 

CL/max 

1.1125 

max.  Lift  coefficient 

experiment  least  squares  fit 

CDmax 

max.  drag  coefficient 

experiment  least  squares  fit 

Table  2:  Model  parameters 


3.4  Rotational  Flow  Model 

In  order  to  have  a  more  complete  model  we  propose  to  include  rotational  flow.  We  used  a  lumped  parameter  model 
that  separates  the  rotational  flow  into  multiple  solid  body  rotations.  Figure  4  shows  the  lengths  and  which 
are  introduced  in  order  to  account  for  differences  in  axial  and  rotational  added  inertias.  The  flow  up  and  downstream 
of  the  control  volume  is  assumed  to  be  dissipative.  With  the  above  assumptions  we  have  to  consider  two  cases  of 
flow  regimes  depending  on  the  direction  of  the  axial  flow. 
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The  model  also  introduces  a  new  form  factor,  K2,  which  makes  the  steady  state  values  of  the  rotational  flow  in  the 
simulations  match  the  experiment.  The  resulting  angle  of  attack  /3  proposed  in  [10]  includes  the  mean  value  of  the 
up  and  dowiistream  rotational  flow  into  the  this  model. 

a  =  atan2(0.7R  (|(a;i  -  W2)  -  i^prop),Vz) 

Axial  force  and  torque  are  computed  in  the  same  way  as  in  the  axial  flow  model  by  replacing  (11)  by  (15). 


(15) 


Figure  3:  Experimental  and  simulated  fiprop  as  a  fimction  of  time  with  variation  of  I  in  the  simulations 
Torque  step  of  +2.15Nm  at  t  =  0.5s,  reversing  to  -2.15iVm  at  t  =  5.5s,  and  to  O.OiVm  at  t  -  10s. 
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Figure  4:  Lmnp  parameters  for  the  rotational  flow  model. 


3.5  Comparison  of  Axial  and  Axial  plus  Rotational  Flow  Models 

Figure  5  shows  the  results  of  the  model  simulations  for  a  2.15JVm  torque  reversal  and  aetu^  experimental  data. 
The  top  plot  of  Figure  5  shows  the  axial  force  versus  time.  Both  models  tend  to  overshoot  significantly  compared 
to  the  experiment;  neither  model  accurately  captures  the  experimentally  observed  transient  response.  In  particular 
the  expanded  rotational  model  performance  is  worse  than  the  axial  fiow  model.  The  bottom  plot  shows  propeller 
rotational  velocities.  It  also  shows  considerable  discrepancies  in  the  transient  behavior  between  model  simulations 

and  actual  experiment.  ^  j  i  ^ 

We  conclude  (i)  that  the  additional  complexity  of  the  rotational  flow  model  alone  does  not  improve  model  pertor- 

mance  and  («)  both  models  are  incapable  of  accurately  reproducing  the  sharp  transient  response  arising  from  step 
torque  inputs  using  experimentally  determined  parameters.  We  note  that  it  is  possible  to  “manually  tune  a  set  of 
parameters  which  will  cause  these  models  to  more  closely  match  the  transient  response  of  a  particular  state,  e.g. 
axial  force,  but  this  will  necessarily  cause  mis-match  in  other  states,  e.g.  flow  velocities. 


4  A  New  Method  for  Generating  Lift  and  Drag  Curves 

Despite  the  addition  of  rotational  fluid  dynamics,  the  thruster  models  described  in  the  previous  sectmn  do  not 
accurately  predict  actual  thruster  transient  response.  This  discrepancy  motivated  us  to  examine  the  validity  o  t  e 
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Figure  5:  Axial  Force  versus  time  (top)  and  propeller  rotational  velocity  versus  time  (bottom)  simulated 
using  rotational  and  axial  model  compared  to  experimental  data.  Torque  step  of  +2.15iVm  at  t  =  0.5s, 
reversing  to  — 2.15iV?7i  at  t  =  5.55,  and  to  O.ONtki  at  i  =  IO5. 

the  sinusoidal  lift/drag  curves,  (10),  [8,  16],  employed  in  the  previous  section’s  models.  Previous  reports,  [8,  16], 
assumed  sinusoidal  lift/drag  curves,  (10),  yet  were  unable  to  directly  verify  these  curves  because  their  experimental 
setups  did  not  include  precision  3-D  fluid  velocity  instrumentation.  We  have  instrumented  our  new  thruster  test 
facility,  [4],  to  precisely  measure  the  variables  (thrust,  torque,  prop  velocity,  and  3-D  fluid  velocity)  necessary  to 
experimentally  determine  lift/drag  curves. 

4.1  Lift  and  Drag  Curves 

Our  thruster  model  assumes  lift  and  drag  are  functions  of  the  square  inflow  velocity  (8)  and  dimensionless  lift  and 
drag  curves.  Figure  6  shows  the  lift  and  drag  coefficients  computed  firom  actual  experimental  data  (at  three  different 
trust  levels),  and  compares  them  to  the  previously  assumed  sinusoidal  lift/drag  curves,  (10).  The  figure  shows  (z) 
a  substantial  discrepancy  between  the  the  experimental  lift/drag  curves  and  the  sinusoidal  lift/drag  curves  and  (^^) 
the  experimentally  determined  lift/drag  curves  do  not  vary  significantly  with  torque  level. 

4.2  Hybrid  Simulation 

The  problem  of  experimentally  generating  accurate  lift/drag  curves  is  complicated  by  two  problems  with  real-time 
measurement  of  thruster  fluid  velocity.  First  is  the  need  to  measure  fluid  flow  velocity  at  the  actuator  disk.  With  our 
acoustic  floppier  flow  instrumentation,  we  were  able  to  measure  fluid  flow  0.1m  up  and  downstream  of  the  thruster  s 
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Figure  6:  Experimental  measured  and  sinusoidal  lift  (upper  plot)  and  drag  (lower  plot)  versus  angle-of- 
attack. 


propeller.  Second,  the  turbulent  flow  ejected  from  the  thruster  has  a  high  variance  in  velocity.  Both  problems  limit 
the  accuracy  of  experimentally  determined  lift  and  drag  curves. 

To  address  this  problem,  we  have  devised  a  novel  “hybrid”  technique  for  generating  lift/drag  cu^es  utilizing 
experimental  data  only  for  (i)  commanded  torque,  (ii)  measured  torque,  (in)  measured  thrust,  and  (iv)  measured 
propeller  rotation  velocity.  First,  these  experimentally  measured  signals  are  used  as  the  input  to  a  simplified  thruster 
model  which  estimates  fluid  velocity  as  follows: 


For  Vz  >  0: 
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For  Vz  <  0: 


Second,  using  the  simulated  fluid  velocity  obtained  from  the  “hybrid”  simulation  and  the  expenmentd  propeller 
rotational  velocity,  we  compute  the  angle  of  incidence  a  and  angle  of  attack  /?  ^  in  (15).  Third,  inve  “S  (  ) 
compute  lifr  and  drag.  Finally,  the  corresponding  lift  and  drag  curves  can  now  be  computed  over  ^ 

of  attack  Figure  7  shows  the  generated  curves  for  diflerent  torque  commands  and  compares  it  to  the  old  hft  and 
drag  coefiicie^s.  The  curves  represent  lift  and  drag  for  one  torque  reversal  by  using  data  from  the  reverse  tr^sition 
we  are  able  to  capture  the  asymmetry  of  the  thruster.  These  curves  serve  as  a  template  for  new  curves  shown  in 
Figure  8.  These  “hybrid”  lift  and  drag  coefiicients.  Figure  8,  can  be  used  in  the  full  dynamic  thruster  simulation  as 

a  lookup  table. 
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Figure  7:  Hybrid  and  sinusoidal  lift  and  drag  coeflficients  versus  angle  of  attack. 
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Figure  8:  Synthesized  lift  and  drag  coeflflcient  versus  angle  of  attack. 


4.3  Comparison 

Figure  9  shows  the  results  of  (a)  actual  experimental  data,  (b)  a  full  dynamic  thruster  model  simulation  (with 
rotational+axial  flow)  using  the  synthesized  lift  and  drag  functions,  and  (c)  correspondmg  siinulations  usmg  the 
sinusoidal  lift/drag  functions.  Figure  9  clearly  shows  the  improvement  of  the  new  /i(/3)  and  /r»(p)  curves,  ihe 
model  using  the  hybrid  lift/drag  curves  predicts  the  experimental  data  with  far  greater  accuracy  than  the  model 
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Figure  9:  Axial  force  versus  time  (top  plots)  and  propeller  rotation  velocity  versus  time  (bottom  plots). 
Plots  show  (a)  experimental  data,  (b)  model  simulation  using  new  lift  and  drag  coefficients,  and  (c)  model 
simulation  using  sinusoidal  coefficients.  Torque  step  of  +2.15iVm  at  t  =  0.5s,  reversing  to  —2.15Nm  at 
t  =  5.5s,  and  to  O.OiNTm  at  i  =  10s. 

employing  sinusoidal  lift/drag  curves.  Moreover,  the  new  lift/drag  curves  also  enable  the  thruster  model  to  capture 
the  forward/reverse  asymmetry  observed  in  the  thruster  experimental  data.  u-  j-  •  i 

Although  the  hybrid  lift /drag  curves  were  computed  using  experimental  data  from  a  single  a  2.15iVm  bi-directional 
step  test.  Figure  10  demonstrates  that  the  resulting  model  is  highly  accurate  for  a  variety  of  torque  levels. 

5  Conclusion  and  Future  Work 

We  conclude  the  following: 

1.  Finite-dimensional  thruster  models  using  either  axial  flow  (Section  3.3)  or  axial-hrotational  flow  (Section  3.4) 
with  sinusoidal  lift/drag  curves  do  not  accurately  reproduce  experimentally  observed  transient  response. 

2.  Experimentally  observed  lift/drag  curves  are  not  simple  sinusoids. 

3.  A  novel  technique  to  experimentally  determine  actual  lift/drag  curves  is  reported. 

4.  Incorporating  experimentally  derived  lift/drag  curves  into  the  axial+rotational  flow  thruster  model  results  in 
highly  accurate  correspondence  between  model  and  experimental  performance  for  both  transient  and  steady- 
state  operation. 
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Figure  10:  Comparison  of  experimental  data  and  model  simulation  using  new  lift/drag  curve  for  different 
thrust  levels.  Propeller  rotation  velocity  (top),  Thrust  (middle),  and  Torque  (bottom)  versus  time. 


A  number  of  questions  remain  unresolved.  The  hybrid  modeling  approach  must  be  tested  and  verified  with  (Merent 
thrusters  and  thruster  configurations.  At  present,  all  model  parameters  are  computed  off-hne  from  experimental 
data  -  an  on-Une  estimation  technique  would  be  useful.  Finally,  a  closed-loop  thrust  control  algorithm  incorporating 
a  form  of  on-line  adaptive  parameter  estimation  with  a  minimum  of  instrumentation  might  enable  improved  model- 
based  thrust  control  in  practical  underwater  vehicle  applications.  Such  a  thrust  controller  may  enable  improved 
closed-loop  positioning  and  tracking  for  marine  vehicles. 
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